The d = 6 trace anomaly from quantum field 
theory four-loop graphs in one dimension. 

Agapitos Hatzinikita^ 

University of Crete, 
Department of Applied Mathematics, 
L. Knosou-Ambelokipi, 71409 Iraklio Crete, 
Greece 
and 

' University of Athens, 

■ Nuclear and Particle Physics Division, 

H ' Panepistimioupoli GR-15771 Athens, Greece 

and Renato Portugal ^ 

Laboratorio Nacional de Computagao Cientifica, 
^ _ Av. Getulio Vargas, 333, 

m ; Petropohs, RJ, Brazil. Cep 25651-070. 

o 

cn 

O I Abstract 

We calculate the integrated trace anomaly for a real spin-0 scalar field in six dimen- 
sions in a torsionless curved space without a boundary. We use a path integral approach 
for a corresponding super symmetric quantum mechanical model. Weyl ordering the corre- 
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O . sponding Hamiltonian in phase space, an extra two-loop counterterm g yR + g^^Tf^-Tij j is 

produced in the action. Applying a recursive method we evaluate the components of the 
metric tensor in Riemann normal coordinates in six dimensions and construct the inter- 
^ ' action Langrangian density by employing the background field method. The calculation 

' of the anomaly is based on the end-point scalar propagator and not on the string in- 

spired center-of-mass propagator which gives incorrect results for the local trace anomaly. 
The manipulation of the Feynman diagrams is partly relied on the factorization of four 
dimensional subdiagrams and partly on a brute force computer algebra program devel- 
oped to serve this specific purpose. The computer program enables one to perform index 
contractions of twelve quantum fields (10395 in the present case) a task which cannot 
be accomplished otherwise. We observe that the contribution of the disconnected dia- 
grams is no longer proportional to the two dimensional trace anomaly (which vanishes 
in four dimensions). The integrated trace anomaly is finally expressed in terms of the 
17 linearly independent scalar monomials constructed out of covariant derivatives and 
Riemann tensors. 
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1 Introduction 



The trace anomaly, namely the breaking of classical conformal invariance of gravity actions 
under Weyl rescaling of the metric: 

gij{x) ^ n{x)gij{x), (1) 

has a long history with numerous applications and implications to high energy physics, general 
relativity and statistical mechanics. The literature on this subject is vast and in this brief review 
we will only concentrate on those aspects which are associated with the present problem. For 
a historical review the reader is advised to consult [|I|. 

Interesting different methods have been developed to investigate and calculate Weyl anomaly 
in four and higher dimensions. The authors in were able to express the integrated trace 
anomaly in four dimensions as a linear combination of two invariants, the square of Weyl tensor: 

a,,ic^^'^ = 4,, - 24 + (2) 

and the only parity even candidate: 

Ci — Kijki ri — K^jf^i — + K [6) 

which is proportional to the well-known Gauss-Bonnet topological density and * denotes the 
dual. To be more concrete, the gravitational contribution to the anomaly depends on only two 
constants (call them a and /3) and is expressed as 

g'^ < T,, >= a (C'^^'Q^M + ^ni?) + (4) 

The numerical values of the constants are: a = ^^ttj, 3 = — ^t^ttt and can also be found 
using the Feynman diagram scheme in [0. 



It was soon realised |^ that these invariants were manifested in the t-independent 62 co- 
efficient of the Schwinger-De Witt asymptotic expansion of the heat kernel of the appropriate 
differential operators. The four (and partially the six) dimensional anomaly was later rederived 
by Bonora et. al. Q who established the connection between Weyl anomalies and cocycles by 
relying on a cohomological method and using the Wess-Zumino consistency condition. Although 
these authors explicitly specified the invariants (see Appendix A. 2 for their classification) they 
did not express the six dimensional trace anomaly in terms of these invariants. From the rep- 
resentation theory point of view Fulling et. al. ||^ were also able to determine the number 
of independent scalar monomials of each order and degree up to twelve in derivatives of the 
metric. The first explicit result was given in the literature by Gilkey ^ and later by Avramidi 
1^, who by making some modifications and innovations to the "heat kernel" or "proper time" 
method, presented a new covariant nonrecursive procedure and found the one-loop effective 
action in the presence of arbitrary background fields in six and eight curved space-time di- 
mensions. Later Deser and Schwimmer re-examined the diferrent origin of the topological 
(type-A) versus local conformal scalar polynomials involving powers of the Weyl tensor and its 
derivatives (type-B) contributions to the anomaly in general dimensions. 

In the present work we follow the supersymmetric quantum mechanical approach first pi- 
oneered by Alvarez-Gaume and Witten ^ and used to compute chiral anomalies. According 
to this method the operators 75, V^, x^, •y^ are represented by operators of a corresponding 
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quantum mechanical model, and by turning these operator expressions into path integrals, one 
finds that anomalies of quantum field theories can be written in terms of Feynman diagrams 
for certain sigma models on the worldline. Bastianelli and van Nieuwenhuizen applied this 
method to trace anomalies |T0|. These authors used mode regularization, a scheme widely used 
at the time. Subsequent work by de Boer et. al |]lT| showed how to use time-slicing and gave 
a completely and unambigous derivation of trace anomalies in terms of path integrals, using as 
input Einstein Hamiltonians. In this article we apply the regularization method of [|lT| to the 
calculation of trace anomalies, following the set up of . 



The outline of the paper is as follows. In the next section we present explicitly one of the 
basic ingredients of the background field method, namely the expansion of the metric tensor 
components in six dimensions in Riemann normal coordinates. The method we use is recursive 
and enables one with the application of (P) and (|lT]) to determine the expansion of the metric 
tensor up to the desired order. 

Section 3 begins with a very rapid introduction to trace anomalies from the one-dimensional 
path integral point of view. We write down the interaction Lagrangian density and the propa- 
gators of the fields involved. 

Section 4 is devoted to the calculation of the perturbative expansion depicting at the same 
time the Feynman diagrams associated to each vertex. 

A vital tool in our journey is a computer algebra algorithm which proved to be very efficient 
especially in finding the contribution of Je vertex with twelve quantum fields. Thus section 5 
is devoted to a brief description of this program. We illustrate it's capabilities by applying it 
to the Ij interaction with eight quantum fields. 

Our conclusions are given in section 6. Several appendices follow to assist to a deeper 
understanding of the technical obstructions the reader might face. 



2 The recursive expansion of the metric components in 
RNC in six space-time dimensions 

Before embarking on the background field method we discuss the expansion of the metric tensor 
in Riemann normal coordinates (RNC). 

RNC have the appealing feature that the geodesies passing through the origin have the 
same form as the equations of straight lines passing through the origin of a Cartesian system of 
coordinates in Euclidean geometry [|12|. Locally no two geodesies through a point V intersect 
at another point, and the power series solution of the geodesic equation is: 

?/' = r^ + E7:y(rU...0^re---e^^ (5) 

k=2 

where {p\-i^i2---i^ ^ the "generalized Christoffel symbols" at the point V and the geodesies 
through V which are straight lines are defined in terms of the arc length s by: 

y' = ^s. (6) 

By induction, one can easily prove that: 

^i^A. ■ ■ ■ di,.,Tl_,,,) = (7) 
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Paraphrasing eq. one can state that all symmetric derivatives of the affine connection 
vanish at the origin in RNC 

In general a covariant second rank tensor field on a manifold can be expanded according to: 



OO 



n=l 



d d 



d 



din dCi2 d^i 



The coefficients of the Taylor expansion are tensors and can be expressed in terms of the 
components i?^„p of the Riemann curvature tensor and the covariant derivatives Dj^Tim and 
^kRmnp- Without much effort one can prove that: 



d{i^di^ ■ ■ ■ di 



«n-l-^ in)k 



/n — 


4) 


\n + 





«n-2-"'j„_iA:i„) 



+ (r 



tn-2a in-lkin) 



I ^ a, a ^ k 
a, a 



(9) 



where the interchange of covariant and contravariant indices act independently and symmetriza- 
tion acts only on i indices. Expression (j^) reproduces for various values of n^ the following 
results: 



2 fci2) 

^ n 

2^{ii^i2ki3) 

31- _^ 2 -I - 



D{iiDi^Di^Rl^i^^^^ D(^i-^R2^i^^^Rl 



1445)0 



-D(j^ ■ • ■ Di, R, 



- [J D{i^Di2Ri^auR-i5i(i)k + ^{ii^i2R'i3ki4^i^i(i)a 

16^ 



The coefficients of(||) can be rewritten as: 
9{iidi^ ■ ■ ■ di^)Tk^k2 = 



^In |13| there is a misprint for the n = 4 case. A minus sign is needed in front of the |-term. 



(10) 
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+ d(iA2 



(11) 



Expressions and (|Tl|) compose the building blocks of the current recursive method which 
produces the following results for different values of n: 



k\k2 



d{hdi2di2,dHdi5)Tk^k2 



kik2 



+ 



+ 



D(ijDi^)Tk^k2 - 2 
D(.i^Di^Di,^)Tk^k2 



^{hkii2)^pk2 + ^1 



9{iA2^i,)kTpk2 + '2d(hTP^^^Di,^)Tpk2 + h 

\ {%ikn2^i3)Tpk2 + h ^ k2 



D(>i^Di^Di.j^Di^)Tk^k2 

^(iAi2^iJ^'ii)ki'^pk2 + '^'^{ii'^i2^i3ki^i4)'^pk2 



+ ^d(^ii^'^2ki 



Ri:,p^^)Tak2 + P^h 



ki ^ k-} 



D(i^ Di^ Di.^ Di^ Di^)Tk^ k2 



10^, 
3" 



~R'{hkii2-^is ' ' ' Di5)^ak2 + ^D(^i-^R'^^i^_^-,^Di^Di^)Tak2 



{ii^i2^i3^iikiisy^oik2 



2 

- -D(^i^R'^^f.^i^Rf^i^)pfak2 + D(^i^R^^^^i^ {^uo^is)^pk2 + a^h 
2 

+ 3 {D{t,Riai3Tpk2 + a^k2) Rt^.^)k^ ±ki^ k2 
D(h ■ ■ ■ Di^)Tk^k2 

+ d{h ■ ■ ■ (^i5^'i(,)kiTak2 + 69(j^ ■ ■ ■ di^T'^^^^^Di^yrak2 

+ 15%, ■ ■ ■ d,,Vl,D,.D,Ak2 + 20%5,,f ^3,^ 



Dif^')T(yk2 



+ 14%rr,,,A3---Ae)T„fe, 

+ 105(, ■ ■ ■ d,,Tl,^ {^^.f'.,^Jpk2 +a^k2) 

+ l0d^.A2^lk, {d,As%aTpk2 +a^k2) 

+ 36%,a,,f-,, {d,Jl^D,,)f,k2 +a^k2) 

+ ^^^^Jt2k, {^^, ■ ■ ■ d.,TP.,)A2 + « - ^2 

+ 24a(,,f-,^ {^,,^^J^aD^,)fpk2 + a^h 

+ i5^^^,Tl,^ {di,Vl^D.,.D,,)fpk2 +a^k2 

+ 15 [%f r^,,a.3f {^^,rl^pf^k2 + P^k2)+a^k2 

+ d(,Slk,D,,---D,,)f^k2 + h^k2 



(12) 



If the second rank tensor with components Tk^k2 is replaced by the metric components gkik2 
then the related covariant derivatives (provided we deal with a torsion free affine connection) 



vanish and the above expressions are simphfied. One could derive for n = 5 the result: 
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On the other hand for n = Q one gets: 



(13) 



k2 



34 



M J^kiiii(i)k2 



ki ^ 



55 - - 16 - - -I 

Thus, plugging into(^ expressions (p!3|) and([T^) we end up with the following expansion of the 
metric tensor in RNC: 



gkik2 



9k\k2 ~l~ 21 2^fcin«2fc2^ ^ 

~l~ ^-^il -^^142*3^2^ ■ ■ ■ 



1 6 

4!5 
1 4 

5!3 
1 10 

11 

2 

+ o(r---r). 



Di^Di^Rk^i^i^k2 + (^Rkiiii2mRi.^ij^k2 



Ai ■ ■ ■ Di^Rkj^i_^i^k2 + 2 \ Di-^Rk^i2i2.pRi^i^k2 + ^1 ^ ^2] 1^*^ ■ ■ ■ 1^ 

17 

Ai ■ ■ ■ Di^Rkj^i^iQk2 + (Ai A2-^A;ii3i4p-Rfgi6fc2 + ^1 ^ ^2 



«5 



+ Ai-Rfcii2i3pA4-Rfgigfc2 tr'R>^ihi2pRi3i4l'RisiQk2 ^ ^ ' ' ' ^ ^ 



(15) 



The expression (P^D will play a crucial role in calculating the contribution steming from the 
quadratic Christoffel symbol term in the interaction Lagrangian. It is also in perfect agreement 
with the closed formula for these coefficients which are encoded in the integral representation 

of ra. 



3 The integrated trace anomaly 

Anomalies in (even) n-dimensional quantum field theories^ are expressed in the Fujikawa [p!5| 
approach as 0: 

Anw = lim Tr {^Je'^'^) (16) 

where JT" is the Jacobian = fix)(p of the fields 0(?/) and the regulator TZ for consistent 

anomalies is uniquely determined [Q. For local Weyl anomalies and for real scalar fields one 

^ For a manifold having odd dimensionality one cannot form a scalar out of an odd number of derivatives. 
■^We perform the usual redefinition of the scalar fields (f> ~ giip which leaves the Jacobian invariant under 
the similarity transformation g 3 . 
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finds that: 

n = H - ^n^^R{x) 

= ^9~H^)Pi9'^i^)9H^)Pj9~H^) ~^^'^CR{x); g{x) = detgij{x). (17) 

The second term in (|T^ is the well-known improvement potential term and the dimensionless 
coefficient ^ = takes the value ^(n = 6) = ^ in the present case. The Hamiltonian 

is Einstein invariant which means that H commutes with the generator G{^) = ^{piC{x) + 
d^)Pi) of the infinitesimal target space diffeomorphisms Xi Xi + C,i{x). 

Consider a spin-0 field which lives on an n-dimensional compact Riemaniann manifold 
equipped with its standard (metric-compatible, torsion-free) connection and having no bound- 
ary. In addition decompose the paths x^{a) into a constant part Xq satisfying the free field 
equations and a quantum fiuctuating q^{(T) one vanishing at the time boundaries. The Weyl 
anomaly for a real spin-0 field may represented by the Euclidean quantum mechanical path 
integral [|1^]: 



Anw{s = 0,n) = lim Tr (/(x)e-^^) 

= n ((2^ / n < >) (18) 



i=l 

1 cint cint < cint „„ J 4. 



where f{x) is an arbitrary function, — ^5"*"* = S"];"* + S^*' and ^: 
Sr = ^ fj9iA^o + q)- 9'.3{xo)] {q'q' + b'c> + a'a^) da 

1 fo n - 1 - h i2 is 

+ ^ (^Di-^Di^Rii.^i^j + —Riij^i2kR i^iij^^ ^ ■ ■ ■ q'^ 

+ ^ • ■ ■ Di^Ru,i,j + 2{Di,Ru,i,kR%,, + i ^ j)) " " " 



— ^Djj • ■ ■ Di^Rii^i^j + — (^Di-^Di^Rii^i^kR^i^j + ^ ^ j) + —- 



+ -Ru,^,kR%.,,R\,.Jq'' ■■■q''+ 0{q^^ ■ ■ ■ q^^] 



(q'q^ + 6V + a'a^) da (19) 



-(3h^-f \^{l-A^)R + g^^T\j:^ 



Ij 



{xq + q)dT 



-Phi 





-1 



[1-40{R + q''D,,R + ^q'^q'Wi,Di,R + ^q'^ ■ ■ ■ q'W,, ■ ■ ■ D,,R 



q''D,,---D,,R]+-g^^d,j:\,d,,T%q'-q 



«2 



^The expectation value < • • • > means that all quantum fields must be contracted using the appropriate 
propagators. A bar over the various geometrical quantities indicates that they depend exclusively on Xq. 
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1 - • • - - 1 

hi-i ^^3^ ik'-'tA^ Ij ik'^ia'-'ii'- Ij 



3! 

+ 0{q''---q 



14 



dr. (20) 

As emphasized in the above expression is the continuum hmit of a rigorous discrete result. 

In the action (|19|) (6*, c^) and a* is a set of anticommuting and commuting Lee- Yang ghosts 
respectively [jlOl- Their existence is imposed by the integration over the momenta p*(r) in the 
transition from phase space to configuration space. A measure factor ^/g is then produced at 
each point of the path and by exponentiating it, introducing the Lee- Yang ghosts, we are led 
to a perfectly regular term in the action. The existence of ghost fields also removes ultraviolet 
divergencies at higher loops and as a consequence all integrals are finite. 

The non-covariant fi^ terms, R and FF, which are essential for the general coordinate in- 
variance of the transition element are created by Weyl ordering the Hamiltonian H of 



(17). The contribution of these terms to the trace anomaly is found by first Taylor expanding 
them in RNC and then substituting the partial derivatives of the Christoffel symbols by the 
polynomials of R and DR: 

{i2k) = -g-^ ii2k)il (21) 

3 



(23) 



The symmetrization of the various indices in the above identities is understood with the in- 
clusion of a ^ factor. All other extra terms in (|I9|) and (^) are produced by expanding the 
metric in Riemann normal coordinates with the help of ([T5|). 

Propagators are derived in closed form through the discretised configuration space path 
integrals via a midpoint rule. In this way ambiguities arising from products of distributions are 
resolved. Taking the continuum limits of the propagators one can read off the Feynman rules. 
One then finds that 6{a — r) is to be considered as a Kronecker delta dij and the propagators 
depend on the discrete Heaviside function On = |. The continous two point Green function 
may also be determined by the equation: 

A(a - r) = S{a - r) (24) 



subjected to the boundary conditions on the interval [—1,0] (end-point approach): 

A(0, r) = A(-l, r) = A(a, 0) = A(a, -1) = 0. (25) 
The Feynman propagator is then formally found to be: 

A{a,T) =a{T + l) + {T-a)e{T-a). (26) 
The propagators of the various fields are proportional to (3h and given by: 

q\a)q^iT) = - ^hg'^ A{a , t) (27) 
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q\o)q^{r 

q\(yW{r 



-Phg'^la + e{T - a)] = -(3hg'^A'{a, r) 
-Phg'^lr + e{a - r)] = -phg'^ 'A (a, r) 
-Phg'^il - 6{a - r)] = -phg'^ 'A* (a, r) 
-2[5ht'5{a - r) 
+(5hg'^{a -t). 



(28) 
(29) 
(30) 
(31) 
(32) 



All other Wick contractions of the fields vanish. Regarding the vertices they may be read 
directly from the continuum Sint given by (|19D and (pO|). 

4 The Feynman diagrams and the associated contribu- 
tions 



Perturbative expansion of < e 



>, keeping only terms that cancel the {13%) ^ factor in the 



measure of the trace anomaly, provides us with the following distinct interactions as well as 
connected and disconnected diagrams having at most four loops |: 



1. 



17 



5«6J 



+ j.Riiii2kR i^i^lR i^iaj 



q'^--- q'^ {cfq^ + &V + a 



[a)da > 



(33) 




1 34 
120 5 



7! V/—"/ 1 120 \ ~ ' "^a. — ^ ' 60 



Daft 1—1 D I 3 pabcdi—i D i o pafe 7~)c r^d p i o TDabcd T\ TD 



1 11 
120 2 



_lQ 7~) TDabcd r-ie p i q r-ia rybcde r-i p 

+Ol>'ait rteftcd + J^eJ^bcda 

p D o. c JD |7p p a6 pe fed 

— J^abcdJ^e f J^bedf + 2J^abcdJ^ef ^ 



DaRbc) + ^DaRbcD'^Rdcab + | (DaRbcde) 



^ 120 5 



I? D'' Z?ic I n Do,b p edc p 



bcde 



^The recipee the diagrams have been drawn is: A hne closing on a single vertex is represented by a Green 
function G(tT, cr) with one integration variable while a line connecting two vertices stands for G{a, r) and indicates 
two integration variables. 
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1 /-o rn 



4! 



(1 - 40 A, • • ■ D,,R - \r\,.,' d,,V\,d,j\ 



1 



+ 



3! y ' 



[a)da > (34) 



QQ = - {(5nf 40 \?>n''R - 2D,R-^D,R - 2R-''D,D,Rj 

ill Dab D cde r> i 1 TDabcd D e f D i 13 Dated D D \ 

~ 540 I 2-"- ^bcde+3^ c ^bedf + ^ab ^cdef J 

1 7~)a Dbc D d I 1 f T-)a pbc\ ^ 5 r-) Dabcd 7~) pe 
"1728^ ^ ^d-^ab c + 17280 V J ~ 2304^^ A^ bed 



^^^D-R'W.R,, - ^ (D'^R^'^-f - ^^D-R^'^'^D,R,,aa 



1259 Dabcd D D i 47 pab p cde p 1 Dabcd D n P 
-"-ab -n-cde/ + 86400 ^hcde — J^d^^bR 



518400 



23 Dabcd /~) 7~) pe 7 Dabcd r—i p 7 Dabcd n n Re 



4800""' ^d^e^ abe 320O 

-Rf'^^R ^ ^ R 



1600' 



129600 



1 fl^ 

2 U^. 



i-bed/ 



rO 



< 



-1 J-1 



127 



X 



+ 



X 



5! 



{<y) 



[r) 



(35) 





+ 



2 144 



^DaRbcdeD^R^'^'^"' + JoDaR°''"^'^D'^Rebcd + ^ {DaRbcde) 

— -^DaRbcDdR'^'^^'^ + ^ {DaRbc) — J§DaR°'''DhR — j^DaR"^'^ D^Rcb — (-Da-S) 

I 1 / 108 p r— 1 Dabcd 1 72 p r— 1 pab 3 pr—i p 1 432 7-1 r-jb Dacde p 1 432 r-ia 7-) pbcde p 
+ 5! I 15" -"'"''"i + 15 -n-ab — 2-KU/£ + -gQ-i^ai^ it Kbcde + J^b^ ^acde 

I 144 7-) 7-) Dbead p i 144 r> D D Dadbc QDD 7~) ffob 1 4 p / p \ 192 p pbcde pa 
+ 75-i>'a-t^b^ -Kcd + -^JJaDbJ^cdi*' — onUaJJbJ^ + 3-K \J^ab J ^^abcdJ^ -Kg 

I 64 p pb d Deafe 224 p jbabef Dcd 64 p pb pea \ n f) ( D \' 

+ 15 -f<abcd-K e /""^ i^^abcd^ ^ ef ~ T^^ab^c-^ + ^^{^^abcd) 
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4. 



h = 



X 



/-o r 1 



12 
(r) 



(1-40 



+ 



X 



(r) 



9 '^is'- mk^iA^ In^ y 



+ ^(1-40 



X 



(r) + -(l-40 



Dj. Di.y Rri 



9 



{a)RdadT > 



(36) 



5. 




[00^0]-00+ 



- -(/^^)^|| l|o^^°-^ ~ llii-^ (Robed) + ilo ^DaR'^'^DbR + ^ (-Da^) 



~600 



■^ROR + ^RDaDjjR"'^ — ^R [Rab) — J^R (Rabcd) 



+ 



q'^q'^Di,Di,R {a)R 



+ 2{i-Ai)Rrdi,v\AJ^\jq''q'' 



Di.Rq'^ {a) D^^Rq 



{T)dadT > 



(37) 



+ 



o 



X • 



{te) {^DaRD"-R + \R^R) - ^^R{Rabcd) 



6. 



^-0 rO 



1296 \j3h^ 



< 



Rnn,jq''q'' (gY + 6V + aV) 



X 



X 



Rkisuiq^'q"' + b^c^ + ol'a^) 



Rmi,ienq''q'' (^^ + b^^c^ + a-a") 



{p)dadTdp > 



(38) 
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oo oax) + [ o X oo ] + cx^ 



-[OOxQOO]- 





+ 




( ^ I ^ P3 I 37 uabpcdD 3_ / p \^ f? 40751 pgb pcde p 

y/Jit-J ]^296 64 360 J^acbd j^g \^ahcdj ^ ' 10080 a bcde 

]_D ( V ^ ^ 20911 TDabcd r)cf p , 409 Dabcd pe / p 137 pab pc p 

8 V ''/ 10080 ab-'^'cde/ "T" 280 « c^bedf 2016 



/7 = 



X 



11 /-O /-O /-O 



T)RdadTdp > 



(39) 



X • 



X • 



JqR^ ~ I {Rabf — \ {RabcdY 



Is = T^Ph (1 - 40' < £ £ £ [i^^ni..-^^^?^^ + b'c^ + 



768 

X RRdadrdp > 



[<y) 



(40) 



X 



OO =m'M^Y\R' 



9. 



1 
3! 



-1371 



40 < J J J ^ RRRdadrdp 



> 



(41) 



1 3 



30^ (4)' 



Some comments are in order 

• There are Feynman diagrams that satisfy the factorization property according to 
which a diagram breaks down into simpler subdiagrams. The four dimensional 
building block diagrams are depicted with their corresponding contributions in what 
follows: 
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Let us examine now the contribution oi terms. In n—2 dimensions the trace 
anomaly is: 

1 



a2 — 



247r 



-R. 



(42) 



The situation changes in four dimensions in which the R^ terms, steming from the 
interactions: 



+ 



X 



+ 



—^Rm.jR < /_° Q^'q'' {'i'q' + + a'a^) da > 



72 {pn) 



R 



[T)dadT > 



1152 



(43) 



and represented by the following diagrams: 

OO --i-m\^:fR^ 
OOO =M^)lK^)'^'-i^(^-)^ 



\2 1 (l. 
2 I 24 



R' 



produce a vanishing result. 

The same observation holds for the R terms created by the interaction: 
Ir = ^\Rii,i,j< j\''q''{q'q^ + h'cP ^ 

- P^^R (44) 
13 



and associated with the Feynman diagrams 



24 24 



One making use of the identities (9)-(24) of the appendix and the 17 hnearly inde- 
pendent terms hsted below: 

it , n[nab) , ri\riabcd) , riabri^n, , it it itacM 

p rjacde r>b p jyabef jycd rjabcd pe / p Z?n D 

itafeit it ^^g, Kabcd-tt -n- e/i -f^^ a c^bedf , rLUK 

Rab^R"''^, Rabcd^R"'^"'', (DaRbc)'^ , R'^^Dj^DcR^ 

D^^R'^D^Rac, {D'^R'^^'f, a^R, {DaRf (45) 



can deduce for the integrated trace anomaly: 

= 0, n) = Irni (t^^ / ft V^/K) A^;^)) (46) 



1=1 



where 



^ ^ 1 p3 I 7 . , 1 p/o 5293 

(/5;i)3 1296000 129600 ^ ''V 4320O ^ "-^"^J 13063680' 



b r>ac 



I 9287 -^^-g^- 159421 - ^acde^b 18413 - ^abefocd 
-r it it itnrbd — ^ah-Tt it „j„ — — iXnbrd it „ f 

16329600 130636800 "'^^ 26127360 

661 - - 3-- 

+ R'^^'^R^ ^Mbedf + ROR RabOR"'' 

362880 " ^ 21600 5600 

191- - 7 - - l/-\2 

:RabcdaR'''"' + ^D-R'WbRac - ^T— t: {D^R'""'' 



483840 8640 20160 

^D'i? - -^R^'DbD^m - (DaRbcY - f/^a/?)' . (47) 

100800 5040 " 120960 V 604880 V ^ ^ ^ 

5 The computer algebra program 



Vertices Ji to Jg were calculated using the Riegeom package |T8[ , which is a Maple package for 
manipulating generic symbolic tensor expressions in the context of Riemannian geometry. In 
this section we show the process of calculating the vertex J7. The lines beginning with ">" in 
typewrite font are the input in a Maple worksheet. A colon at the end of the command hides 
the result. We start loading the package 
> with (Riegeom) ; 

defining Christoffcl(Ganima), Riemann(R), Weyl(C), Ricci(R), LeviCivita(eta), 
TraceFrecRicci(S) for Dimension — 4, CoordinateName — X, MetricName = g 

[absorbg, changedumind, cleartensor, codiff , coordinate, definetensor, dimension, 
expandcodiff , Iptensor, Itensor, metric, normalform, off' , on, printtensor, replace, 
simpLC , simptensor , sreplace, switches, symmetrize, symmetry, tdijf] 
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Next command reads file "Vertex.mpl" that contains procedures written specifically for the 
calculation of vertices Ii to Ig using Maple programming language and Riegeom commands. 

> read( 'Vertex.mpl ' ) : 

We setup spacetime dimension using Riegeom command dimension. 

> dimension(6) : 

Next command enters the tensor coefficient of vertex J7. 

> tensor_coeff := printtensor(R[-mul, -il, -12, -nul] *R[-mu2, -13,-14, -nu2] ) ; 

tensor _coejf := R^ims v\ i3 n u2 

Prlnttensor is the Riegeom interface command. Indices with minus sign are covariant and 
with plus sign are contravariant. In the next command we enter the field terms. We use a 
Maple list (which preserves order) instead of an expression written in terms of the commuting 
product operator. In this form we have full control over the order of the fields. 

> L := prlnttensor ([ [q[ll] ,q[12] ] (slgma) ,[ [qdot [mul] ,qdot [nul] ] , 

> [b[mul] ,c[nul]] , [a[mul] ,a[nul]]] (slgma) , [q[13] ,q[14]] (tau) , 

> [ [qdot [mu2] , qdot [nu2] ] , [b [mu2] , c [nu2] ] , [a [mu2] , a [nu2] ] ] (tau) ] ) ; 

L:= [[q ^\a\q^'{a% 

[[qdot '^i(a), qdot ''\a)l [h ^\a),c ^\a)l [a ^\a),a (a)]], 
^'(r),g^^(r)], 

[[qdot ^\t), qdot ^\t)1 [h ^\t),c^\t)1 [a '^^(r), a '^^(r)]]] 

Next command finds all independent Q index configurations that contribute to the final result. 

> Lie := lnd_conf lg( [11 , 12 ,mul ,nul , 13 , 14,mu2,nu2] ) : 

> N : = nops (Lie) ; 

:= 105 

In the case of 8 field indices, there are 105 independent index configurations. We show the first 
3 ones. 

> for 1 to 3 do evaln(Llc [1] )=Llc [1] od; 

Lic\ = [il , i2, fil, ul, i3, i4, /i2, z/2] 

Lic2 = [il , i2, /il, ul, i3, /i2, i4 , v2\ 

Lics = [il , i2, /il, ul, i3, z/2, /x2, i4] 

Next command is a loop over the 105 index configurations. The results are stored in a table of 
results called res. 

> for 1 to N do 

> expr := WlckContractlons(L, Llc[l]): 

> res[l] := Vertex(tensor_coef f *expr) : 

> od: 

Command WlckContractlons performs the Wick contractions given by eq. (27) to (32). We 
^The number of independent index configurations is given by where (2-)! represents the possible 

permutations of index pairs and 2^ the permutation of indices within the pairs. 
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show an example with the fourth index configuration. 

> expr := f actor (WickContractions(L, Lic[4])); 

expr := hhar^ g'^ g ^2 gf^^i3 ^ ui ^4 ^ (1 + ^) + Heaviside(a - T)f 

Command Vertex simply multiplies the result of Wick contractions to the tensor coefficient, 
simplifies the tensor expression, and isolates the terms to be integrated, since we noticed that, 
for the most complicate vertices. Maple spends more time integrating Dirac delta and step 
functions than simplifying the tensorial terms. 

> Vertex(tensor_coef f *expr) ; 

[/?^ hbar^R"^ ''^ R^2 ui, cr {I + a) {t + Heaviside(cT - r))^] 

Next command adds the results of all index configurations after integrating over variables a 
and r. Functions int_tau and int_sigma replace the usual Maple integrator, since Maple int 
command fails to return the correct value for complicated vertices. 

> f inal_value := simptensor (add(res [i] [1] * 

> int_tau(int_sigma(res [i] [2] )) ,i=l. .N)) ; 

final _value ■= — hbar^ - - hbar^ R"^ R ^1 .2 - 
16 6 

1 4 4 21 

- j3 hhar R" - ~^ Ru2ii _bo -co 

Collecting ^^Z?"^, we obtain the final form. 

> collect (final_value , [beta,hbar] ) ; 

(77: R^ — - R'^^ R ui u,2 ~ R'^'^ "'^^ Ru2 ii _bo -co) hbar'^ /J"^ 
16 6 4 

For vertices I2 and Iq, the method described above has some extra complications. Vertex Jg 

has a large number of index configurations (A^ := 10395) and vertex I2 uses huge side identities 

obtained in Riemann normal coordinates as described in appendix A.3. 

This computer algorithm can easily be extended to higher dimensions to predict the trace 

anomalies in 8 and 10 dimensions. Of course there are strong limitations implied by the 

exponential growth of the problem which can be relaxed by increasing the available computing 

power. The tensorial upper bound of the program is reached by the product of eight Riemann 

curvature tensors but such a case is beyond the scope of the present work. 

6 Conclusions 

In this article we calculate the integrated trace anomaly for a real spin-0 scalar field living on a 
curved six dimensional manifold. This is achieved by relying on a recursive computation of the 
metric tensor components in Riemann normal coordinates. One can use the general formulae 
d^) and (p!TD to reach the desired order of metric expansion induced by the dimensionality of 
the manifold. Adopting the path integral formalism of quantun mechanical non-linear sigma 
models, we evaluate all the vertices and the corresponding Feynman diagrams that contribute 
to the present anomaly. A computer based program is used to perform the otherwise inevitable 
task of integration over distributions and contractions of the various tensors involved. The final 
result derived in this process involves 17 scalar monomials consisting of covariant derivatives 
and/or Riemann tensors. 
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Appendix 



A.l Useful identities 



We consider a Riemannian (or pseudo-Riemannian) manifol equipped with its standard (metric 
compatible, torsion free) connection. The Riemann curvature is defined as: 

R\ca-dj:\a + ^\rec-c--d (1) 

possessing the famihar symmetries: 

• antisymmety 

Rabcd — —Rbacd — —Rabdc (2) 



pair symmetry 



cychc symmetry 



Rabcd — Rbadc 



Rabcd + Radbc + Racdb — 

Bianchi symmetry and the related identities 

Rabcd;e ~l~ Rabec,d ~l~ Rabde,c — 

Rbc;a ~ Rab;c + R"^ bca;d ~ ^ 
1 



The Ricci curvature and scalar are: 

Rnh — R'^ 



Rb - i^KR 



R^R^ 







t-ab aqbi 

In this paper the following identities have been exploited to simplify our expressions: 

D T>acbd 
^abcd^ 



-ry ryecbd 
J^abcd^ 



Ra{bc)dR^''"'^''^ Rf% 



2 abed 




D jyebcd 
-J^abcdJ^ 


lo 

Q J^abcd 
O 




-n^R 
2 




-n^R- 
2 







R\iR\I + ^R^'^^Ref' 



(3) 



(4) 



(5) 
(6) 

(7) 



(8) 

(9) 
(10) 

(11) 

(12) 



if -2 

yab Dcd ; 



ia rib 1 



i'acbd 



+ 3 {DaRbcY + T^RabcdPR"''^'^ — '^RabcdR^ af^R 



b T^ac 



ADaRbcD'R 



2 
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di 



I O Dd'b TDcde r> rtahcd d D _l I? r~l D^h (^'^\ 

+ ZK K a^bcde — ^abef-n-cd + ^ab^^ K^'^) 

DaD.D'^D.R'"' = ld'R--^{DaRf -]^RabD''D''R (14) 
D^D.D.D'R-'^ = ^n^R-^(DaRf-2Ra^D''D''R 

- ^DaRbcD^R"^ + Rab^R'''' + '2{DaRbcf (15) 

DaDbD^D'^R"' = ^n'R^^{DaRf -2RahD''D'R 

+ Zriabiic^ — ZK K Kacbd ~ oJJa^bcJ--' -K 
+ Rab'^R"^ + 2 (DaRbc)'^ + -Rabcd^R"'^'^ 

Tjabcd pe c p6e i pab pcde p 1 ryabcd p P 

-it K K + K K a^bcde - ^abef^cd 

D'^R^'^DaRj, = {D.Rbf - D^Rb.D'R"- (16) 

rtabcd T> ryohcd n p i P pacde p6 p pe c pb/ / 1 7^ 

it LfdUbriac — ^-n- ^^abcd + -^^ab^ ^ cde ~ ^abcd^ af ^ de \^ ' ) 

lyabr^cr^ td d D'^^'n Z? Z?'*'' 7") D H -L. Z?'*'' Z?"^ Z? Z?"*^ Z?'^'' Z? M S^ 

-ft U UdrC^ fj — K i-irCab — -K L>aJ-'b^'ri^ J^a^cb ~ ^ -K -Kacbd l-LoJ 

R''''^DdDeR\,, = R^'^^DaDbRac (19) 

DaR'""'''DeR\,a = 2 [(L>„i?6e)' - i^ai^bci^'i?"'^] (20) 

R-''^D,DaR\,, = ^i?"'^"ni?„bcd (21) 

D°^R^^'^DeRabcd — 2 i-^aRbcdef (22) 

DaRD^R = ^UR'^-RUR (23) 
R'^'DbDaR = 2R-'D'DbRac-2R'''R'''Racbd + 2R'''RlRbc (24) 
Proof of (9) By cyclic symmetry, 

p TDacbd p / ryadcb i ryabdc\ p ryadcb , p2 ^'Oc;^ 

^abcd^ — —^abcd + ^ ) — —^abcd^ + ^abcd' l^'^J 

Rename the indices in the first term: 

p ryadcb p r>acdb p ryacbd 

n-abcdJ^ — iiabdc^ — ^abcd^ [ZO) 

yabcd ne ; 



Solve for the desired object. In the same vein we can prove (10) as well as DaR D^R^cbd 

' ^ ,^ „i Tiabcd I — i Ttaticfl. i — i 7 

^abcd 



\ [DaR"'""^) and R^^'^'^UR^^i^a = |i?"^"'^ni?a 

Proof of (11) Expanding out the products, the resulting sum can be expressed in terms of 
the invariants: 

7" p pabe/ pc d 

J^l — ^abcd^ ^ e f 

T p pabe/ p cd 

J-'2 — iiabcd^ ^ef 

T p Z?e c p/ b 

-t^3 — i^abcdi^ fa ^ de 

Li = RabcdR^afR^ed (27) 
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and with the assistance of the identities: 

L, = 
Li, = -L2 

L4 = —L3 + RahcdR^afL^ fd (28) 

one can recover the proposed formula. 

Proof of (12) - (24) With the assistance of ('Tj - (^5^ it is a straightforward excercise to show 
their vahdity. 

A. 2 Six dimensional invariants 

In d-dimensions the polynomials C°'^^'^CcdefC^^ab ^'^d CabcdC'^'^'^'^C^'e/ (type-B anomaly according 
to the geometric classification of 0) are written as: 

r™., f^abcd^ f^ef Tjabcd ije/ p pafece pd , ^ D d2 

InVi - Ly L-cdef^ ^f, - K Kcdef-tt ab - ^ _ 2 ^a.bcd^ -^g + _ ^-^ _ J^^abcd 

8(2rf-3) 3 24(2rf-3) , 

(d- l)2((i-2)3 {d - l){d - 2)^ 

16((j — 1) fj pc , 24 jyabcd p p /r)QN 

(^d _ 2)3 -"-fec-Ka + _ 2)2 -"^ J^ac^bd- [^y) 



RR?„h 



3 pp2 , 2(3d-4) , , 

2(d-l)(d-2) (rf-2)3^"'^ ^'^ 



3r/ 3 

I p pac pbd p pafece pd 
l^^abcdn. K — -riabcdJ^ -Kg 

+ RabcdR^'^^R^/. (30) 

Expressions(p9D and(^0]) reproduce in = 6 the unique already known Weyl invariant polyno- 
mials of 0] constructed only out of Weyl tensors. In higher dimensions the independent ways 
to contract indices among a number of Weyl tensors increases. 

Another invariant made out of covariant derivatives of the Weyl tensor is: 

InV^ = —R^ + SRR°^Rab + 2 RR"''^'^'^ Rated — ^^R°^ RlRca — ^^RabR"^'^ Red 

+ ]^RnR-bR''''uRat + ^R-'>^''nRa,,d. (31) 

The final nontrivial ^ invariant (type-A) we would like to consider stems from the Euler 
form which exists in any even dimension d = 2n: 



E2n = 7 1 r/ i?"^"' A- ■■ A i?"2n-ia2n 

(47r)" n\ Jm ' ' 



^The terminology trivial invariants is adopted to justify the existence of covariant total derivatives of poly- 
nomials over Riemann tensor and it's covariant derivatives. These invariants coincide with variations of all 
independent local counterterms to effective action. 
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where 



1 1 r 



1 



(47r)" n! 



^2nC?V^ (32) 



'-'^n 2n "■^"■2---a2n'^ ^ bite'''-"" 62n-lf'2n 

= ^det(5^0/?--,,,,---i?"^"--\_,,„; ^ = l,---,2n (33) 

is the Euler number which is a total divergence, vanishes in all lower integer dimensions and 
dV = A • ■ ■ Ae^" = ^/—q(P"'x is the volume form. In six dimenions the Euler number becomes: 



InV4 = £q = — 12RR°^Rab + IGR'^^RlRac + '^^R"'^ R^'^ Racbd + 'iRR'^'^'^'^Rabcd 

I A ryabcd JD r> OA Ty<^^ T? ^de p o Tyabcd D e / p 



A.3 The vertex /2 

The key idea for calculating the contribution of this vertex is to express each interaction term 
in terms of polynomials of the Riemann curvature and its covariant derivatives. This can be 
achieved by making use of the identities (|2^) and the symmetries implied by the structure of 
each term separately. 

• With the help of the identity daT^ = daT'^^^) = -lRfbc)a 

one may rewrite the second term 

of I2 as 0: 

• Expanding out the identity d(^adh^c)d — ~\^(a-^b\d\c) S^^- 



Notice that the third term of I2 is symmetric under the interchange of the lower indices of 
the Christoffel symbols namely ^ k and I ^ v. So in(^) we interchange these indices 
and add the two expressions together. We then obtain: 

2 (5n5.2rL + 5n%r[).2 + ^^2%r[)n) 

~~^i'^il^kii2 "I" ^h^kiii "I" ^ii^iki2 ~^ ^i2^ikii 
+ DkR'^2^^, + D^R^ + Ai?!,fc,, + Ai?U,J. (37) 

^Identical results arise if one makes use of the relation 9(af — ~\Rfbca) which communicates with the one 
employed in the text by the relation 9(0? — 0. 
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We would like now to express the second and third terms of the lefthand side of (|37|) in 
terms of di-^di^T''^.^ and covariant derivatives of the Riemann curvature. Starting from: 

^n%f|,)j2 = •^ji^jjf (j^) — (38) 

one has: 



Substituting(^) into(|37|) one gets: 



+ DkR\^ii^ + DkR\^ii.^ + DiRl^^i^ + DiR^^^iJ 

+ 2 ('^^i-^^^zfe + -^n-^Ui + ^«2^Lifc + Aa^Lii) • (40) 

Finally the expression we are going to use for the contribution of the third term in the 
action I2 is: 

f^ii<9i2r'^j = — — {Di^R''i^--^ + Di^R''f^--^ + Di^Rlf^-^ + Di^Rl'^^^^ 
+ DkR\^ii^ + DkR\iii^ + DiRl^f^i^ + DiR\^f,^J 

+ Q{DhRu2k + DiiRii2i + Di2R[i^f^ + Di^R^i^^Y ^^^^ 
and a similar relation holds for the term di.^di^Ti-. 
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